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SUMMARY 

A set of higher-order boundary-layer equations are derived that 
are valid for three-dimensional compressible flows. The equations 
are written in a generalized curvilinear coordinate system, in which 
the surface coordinates are nonorthogonal, whereas the third axis is 
restricted to be normal to the surface. Also, higher-order viscous 
terms which are retained depend on the surface curvature of the 
body. Thus, the equations are suitable for the calculation of the 
boundary layer about arbitrary vehicles. As a starting point, the 
Navier-Stokes equations are derived in a tensorian notation. Then by 
means of an order-of-magnitude analysis, the boundary-layer 
equations are developed. To provide an interface between the 
analytical partial-differentiation notation and the compact tensor 
notation, a brief review of the most essential theorems of the 
tensor analysis related to the equations of fluid dynamics is given. 
Many useful quantities, such as the contravariant and the covariant 
metrics and the physical velocity components, are written in both 
notations. 


1. INTRODUCTION 


For the numerical prediction of the viscous flow about three- 
dimensional aerodynamic configurations, two distinct approaches 
are being developed. In the first approach, the Reynolds-averaged 
Navier-Stokes equations are solved for the entire flow field. In the 
second approach, an inviscid-flow method is used for the prediction 
of the outer flow, whereas near the walls the boundary-layer 
equations are solved. The two solutions are then connected by means 
of a viscous-inviscid interaction technique. 

While the Navier-Stokes approach is more general, the viscous- 
inviscid interaction approach may be more computationally 
efficient. In addition, Navier-Stokes algorithms have been developed 
which incorporate boundary-layer solvers in order to improve 
accuracy near the walls. This is called the Fortified Navier-Stokes 
approach(ref. 1). 


A simple and computationally efficient algorithm for solving the 
unsteady three-dimensional boundary-layer equations has been 
developed recently by Van Dalsem and Steger(ref. 2). In the related 
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boundary-layer equations, the surface curvature was neglected. As 
this poses restrictions in the use of this algorithm for those bodies 
with large curvature, it is necessary to develop a more general set 
of equations which includes the curvature effect. 

In this work a set of higher-order boundary-layer equations are 
derived for a nearly-general curvilinear coordinate system. In this 
system, the surface coordinates can be nonorthogonal, whereas the 
third axis is restricted to be normal to the surface. Components of 
both the covariant and of the contravariant vectors appear in the 
equations because of the nonorthogonality of the surface 
coordinates. To obtain a compact set of equations, the theory of 
tensor analysis is used to derive the equations. The surface 
coordinates will be denoted by (x 1 , x 2 ) = (£ , rj), and the normal 
coordinate by x 3 - £ . 

These equations should be easy to use for the derivation of a 
fast algorithm, which will be used in the Fortified Navier-Stokes 
approach. 

Similar derivations of the boundary-layer equations in curvilinear 
coordinates have been described in the literature, but they are 
generally first-order in approximation. The majority of these are 
based on the classical equations of Squire(ref.3), which are valid for 
an incompressible fluid. Robert(ref.4) has treated the case of the 
compressible fluids in a higher approximation, that accounts for the 
effect of the curvature of the wall. The equations he derived are 
characterized by the use of the "shifters" technique that he 
borrowed from elasticity theory. Thus the various geometrical 
quantities above the wall are expressed as functions of their 
projected values on the wall and the vertical distance. Hirschell & 
Kordulla(ref.5) provide a comprehensive review of this topic and of 
the various methods that are used for the inclusion of curvature 
effects. 

As a starting point, the Navier-Stokes equations wil be derived 
in a fully nonorthogonal coordinate system using the compact tensor 
notation. Then the boundary-layer equations will be developed, by 
means of an order-of-magnitude analysis. For brevity, most of the 
essentials of the tensor analysis and of the theory of 
transformations in curvilinear coordinates are reviewed without 



proofs in the following section. For a thorough review, see the books 
of Aris(ref.6), of Hinchey (ref .7), and of Borisenko & Tarapov(ref.8) . 
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2. BRIEF REVIEW OF TENSOR ANALYSIS 

Throughout this work it will be assumed that the geometric and 
physical quantities of the flow are referenced to a Cartesian 
coordinate system y'(x,y,z). The curvilinear coordinate system into 
which the flow field will be transformed will be designated by 
Initially, the system will be assumed to be fully non- 
orthogonal. The various expressions will be given both in compact 
tensorian form as well as in the analytical form (using the letters x, 

y, z, t|, 0- 


Many of the advantages of the tensor method derive from the 
simplifying nature of the tensor notation and of the summation 
convention. For example, if the product, c, of two matrices a and b is 
considered 


C U C 12 C 13 


C 21 C 22 C 23 


C 31 C 32 C 33 



a il a i2 a i3 


a 2l a 22 a 23 
^1 *31 a 33 



b ll b 12 b 13 
b 21 b 22 b 23 
b 31 b 32 b 33. 


( 2 . 1 ) 


then in the tensorian summation convention the elements of the 
matrix c can be calculated by 


c.. = a, b. . = a., b. . + a.- b«. + a., b,. 
ij ik kj ll lj i2 2j i3 3j 


( 2 . 2 ) 


In the example of equation (2.1), for the definition of the 
elements of the 3x3 matrices a, b, c, it was necessary to use two 
indexes. These entities are called tensors of second order. A typical 
example is the stress tensor. A vector is a quantity specified by 
three numbers, namely its components with respect to a given base, 
whereas a scalar is a quantity whose specification requires just one 
number, in any coordinate system. Scalars and vectors are 
considered as special cases of the object called a tensor of order n. 
In this case the order is zero for a scalar (no index is used) and one 
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for a vector ( one index is required). Commonly-used tensors of 
higher order are the permutation index ej jk , and the curvature 
tensor, Ftp . 

In addition to the number of indexes, distinction must be made 
between indexes appearing as superscripts, e.g., u' , or subscripts, 
e.g., gjj . Superscripted indexes are called contravariant tensors. 
Subscripted indexes, covariant tensors. The case of a mixed tensor 
appears quite often in applications, a typical example being the 
Kronecker delta : 5'j . 

A complete definition of a tensor includes the transformation 
law of its components. Thus, an entity specified by 3 m + n 
components, 

jl...jn 

is called a tensor of order m+n if by transformation from coordinate 
0(yT y 2 > y 3 ) to 0(x 1 * x 2 - x 3 ) > it becomes 


A(x) 


pi. ..pm 


9^ 3x^3/ 3/ ajm 

ay 11 Sy™ dx^ '"dx* y jl - jn 


(2.3) 


More precisely, it may be called an absolute tensor of 
contravariant order m and covariant order n . 

The partial derivatives 3x'/3yi and their inversions must be 
known for the calculation of the elements of a tensor in another 
coordinate system. The transformation operations are easy if the 
metric coefficients are used (see the next section). 


2.1 Bases and Metric Coefficients 


A curvilinear coordinate system xi ( r\, t, ) can be defined in 

relation to the basic Cartesian system y'(x, y, z) 


£ = fl (x,y,z), ri = f 2 (x,y,z), £ = f 3 (x,y,z) 
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Then, through each point in a given region of space there will pass 
the surfaces £ = constant, r\ = constant, £ = constant, which are 
known as the coordinate surfaces. Any two of these constant 
surfaces intersect in a space curve. The three curves that pass 
through a point are known as the coordinate lines of this point. If r 
is the radius vector from an arbitrary origin to a point, the change in 
r caused by infinitesimal displacements along the three coordinate 
curves is 


dr dr dr 

dr = — + — dT| + — d£ 

^ dri d C 


( 2 . 4 ) 


The set of vectors 

dr dr dr 

1 d^ 2 dn 3 dc 

are known as the covariant base vectors of the curvilinear 
coordinate system. The vector ei is tangent to the coordinate curve 
(£) and points in the direction of increasing The base e-j, e2, e3 is 
said to be local, since it generally varies from point to point. In 
general, these vectors are neither perpendicular nor of unit length. 
Analytically, they are defined by the relations 


ej -i +jy^ + k z^ 

e- =i x + jy + k z 
2 T1 T1 11 

e 3 - j y ; + k (2 , 5) 

where (i, j, k) are the unit vectors of the Cartesian system. 

Since the set of the covariant base vectors are noncoplanar, they 
define a parallelepiped whose volume is given by 

V =e 1 .(e 2 xe 3 ) =62.(63X6^ =63^x62) 

A triplet of vectors known as the contravariant base vectors is 
defined by the relations 
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i 

e = 


e 2 xe 3 


2 

e 


e 3 xe i 


3 

e 


e i Xe 2 


( 2 . 6 ) 


It is clear from equation (2.6), that e 1 is normal to the plane 
defined by the covariant base vectors (e 2 , e 3 ), etc. Analytically the 
contravariant base vectors are defined by 

«' =■ ^ + j 'V,. +k 
e 2 =i n. +j H y +k n z 

e 3 => t x +j C, +k C 2 (2.7) 

An increment of arc length along a space curve is given by 

( ds ) 2 = ( dr ) 2 = dr .dr = ej .e^ dx 1 cbd (2.8) 

The nine dot-products e j.ej form a symmetric tensor called the 
covariant metric tensor gy = gjj = ej .ej . Thus the length element is 
given by 

( ds ) 2 = g ;j dx 1 dx j (2.8a) 

From the classical theory of partial derivatives, it is known that 
a volume element is transformed by the equation 

d^ dr| dri = — dx dy dz = J dx dy dz (2.9) 

8(x,y,z) 

where J is the Jacobian. 

In several textbooks on tensors, the volume element alternatively 
is given by 

dV = dx dy dz = yj d^ drj d£ (2. 10) 


where g is the determinant of the metric tensor: 
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g ll 

gl2 

§13 

&21 

§22 

§23 

§31 

§32 

§33 


( 2 . 11 ) 


Thus, the Jacobian of the transformation is related to the 
quantity g by the relation , 

J =_L (2.12) 

V§ 


The contravariant metric tensor is defined in a matter similar to 
the covariant: g'i =gi’' = e'.eJ . Its elements can be estimated by using 
equation (2.7). Alternatively, they are connected to the elements of 
the covariant metric tensor with the relations 


8 = y = g< 8 *, 8 kn- 8 j. 8 km > 
(i, j, k) cyclic (1, m, n) cyclic 


(2.13) 


where G' 1 is the cofactor of g^ in the determinant g . 

Other quantities that are derived from the metrics are the 
Christoffel symbols of the first kind 


r • • , , 1 , 8g ik 8g & d g jJ v 

IlJ ' k,= 2 ( ^ + 3^3x k) 
and of the second kind 

(i'j )-{/,) =r i ‘ i = g B [ij.k], or = g M ■ ~Tj 
,J 11 1 3x k 


(2.14) 


(2.15) 


The Christoffel symbols may be interpreted in terms of the 
variation of the base vectors with respect to the coordinates. They 
appear in the expressions of the covariant derivatives (§ 2.2). 
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In orthogonal curvilinear systems, such as the cylindrical or the 
spherical system, many of the metrics are equal to zero. More 
specifically, in the case of an orthogonal curvilinear system: 


As this relation is simple, the following symbols are usually 
chosen for the nonzero gy : 

N = -J& n » ^2 = J §22 ’ ^3 = J &33 

Finally, in a Cartesian coordinate system (rectangular) the 
covariant and the contravariant base are identical. 

In the next section the usefulness of the various symbols that 
have been reviewed here will be demonstated. The metrics and the 
Christoffel symbols of the second kind are given analytically in the 
appendix for the coordinate system x' {%, t|, Q . 


2.2 Derivatives Based on Tensor Algebra 


If A ik and B ik are the components of two second-order tensors, 
then the numbers C ik , defined by 

C ik = A ik + B ik (2.16) 

are the components of a second-order tensor , called the sum of the 
tensors with components A ik and B ik . Addition of any number of 
tensors of arbitrary order is defined similarly. Tensors of different 
orders cannot be added. The tensors must have not only the same 
order but also the same structure, i.e., the same numbers of 
covariant and contravariant indices in the same places. Subtraction 
of tensors is defined similarly to the addition. 

On the contrary, tensors of arbitrary order and structure can be 
multiplied. The general rule , for multiplication(outer product), is 
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C ijk =A ij B k (2.17) 

mm v 7 

The operation of summing a tensor of order n (n>=0) over two of 
its indices is called contraction. For example, a possible 
contraction for the tensor A^'i is the operation 

Ajf = Aj 1 +A‘ 2 2 +A‘ 3 (i =1,2,3) 

The resulting tensor actually is a contravariant vector. 
Contraction can be performed only on pairs of indices in different 
positions, i.e. , one contracted index must be covariant and the other 
contravariant, in the case of generalized coordinates. Otherwise, 
the result of the contraction will not be a tensor. 

The result of multiplying two or more tensors and then 
contracting the product with respect to indices that belong to 
different factors is called an rnner product. For example, the 
relation (2.17) if contracted, is written 

C ik = A ij B k (2.18) 

Another useful class of operations is the raising or the lowering 
of the index. If Aj is a covariant vector, the contravariant vector 

Ai— g>iAj is called its associated vector. This operation is known as 
raising the index. Its inverse is the lowering of the index by an inner 
product with gy : Bj= g jj B ' . The operation is also applied to tensors 
of higher order. 

It has been mentioned previously (§2.1) that in a generalized 
system of coordinates the base vectors are local, i.e., they are 
functions of the coordinates x 1 , x 2 ,x 3 . Then it follows that the 
differential of a vector includes not only a term that expresses the 
change of the components of the vector itself, but also a term 
required because the base of the coordinate system also varies from 
point to point. The name covariant differentiation is applied to this 
operation. The covariant derivative is denoted by A,j (suffixes have 
been suppressed in the vecior A). The covariant derivative is a 
tensor. Standard formulas are included in some of the reference 


books for the calculation of the covariant derivatives of tensors of 
various orders. Here we are interested in formulas up to the second 
order 
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u. = 

l.k 

5u . 

' . } u . 

i k J 

(2.19a) 

u !k = 

au 1 

^ + 

{ \}u j 

J k J 

(2.19b) 

r^ik 
,1 ~ 

3T ik 

+ 

ax 1 

{ ' , ) T mlt + { k , } T im 

1 m 1 J 1 m 1 J 

(2.19c) 

T 

ik,l 

1) 

. { m ) t - f m } T 

1 i 1 1 1 mk >k 1' im 

(2.19d) 


For the derivation of the flow equations, Ricci’s theorem is very 
useful; it states that the covariant derivatives of the metric tensors 
vanish, i. e., the components of the metric tensors and the 
determinant g can be regarded as constants under covariant 
differentiation. Also, in the case of a zero-order tensor (scalar) , 
the covariant derivative reduces to the partial derivative with 
respect to the coordinates 

rlf 

f i=- (2.20) 

° dx J 

Thus, the covariant derivative of a scalar f is a covariant vector 
with components equal to the covariant components of grad f . 

The divergence of a vector A = A' is defined as the contraction of 
the covariant derivative of A' 

3A* 

divA = A . = + { .* . } A J (2.21) 

" dx* 1J 

The contracted Christoffel symbol may be expressed, after some 
manipulations, in terms of the determinant of the metric tensor 
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i 1 d i — 1 

{■■}=-= — (Jg) ( 2 . 22 ) 

J Vg 0X* 

Substituting the above expression in equation (2.21) , the 
following results are obtained for a contravariant or a covariant 
vector respectively 

div A = -p^-r( A‘7g ) = i_^_(g lJ A 7g ) (2.23) 

Vgdx* Viax 1 J 

Equation (2.22) may also be used for the estimation of the 
covariant derivative of a second-order tensor (equation 2.19c), for 
the case in which the operation of contraction has been applied: 

T !k = ^-^r < V? T " n) + ij k ) T " k (2-24) 

This relation will be used for the calculation of the viscous 
forces in the momentum equations. 


2.3 Flow Quantities Expressed in a Generalized Curvilinear System 


In order to relate the velocity components of the flow in a 
generalized coordinate system to the corresponding velocity 
components in the reference Cartesian system, the relation between 
the differentials of the coordinates must be found. This is easily 
done by application of the chain rule 

dx 1 = dy j (2.25) 

V 

Then, by taking the time derivative of the contravariant vector 
dx' , the following expression is found 
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u i_d^_ax/dx^_3x^ u j (2.26) 

~ ~ 3 y j dt " 3 y j U 

Thus, the contravariant velocity components, U' , are connected 
to the Cartesian velocity components, u', by the components of the 
contravariant base vectors. Analytically, the following well-known 
equations are valid 

U = u £, +v£ + w £ 

V = UTl x +VTl y +WT) z 

W =uC x +vC y +w C z { 221 ) 

Similarly, the acceleration and all higher derivatives are 
contravariant vectors. If the derivative of a scalar quantity is 
considered, the resulting vector is a covariant one. Indeed, by 
application of the chain-rule 


£- = (2.28) 

9x* 3x* 

It is seen that the components of the covariant base vectors 
connect the two flow fields. The scalar f may be the pressure, the 
density or other variables. As an example, the ^-derivative of the 
pressure, p, is given by 

= V* + Vy + V* 

The three Cartesian coordinates have all the physical dimension 
of length, but in general this does not happen in a curvilinear 
system. For example, in the standard cylindrical polars, two of the 
coordinates have the dimension of length, but the third, being an 
angle, has no dimensions. Thus, the contravariant velocity 
components will not all have the same physical components of 
velocity. The same problem appears in all the vectors which are 
related to the flow, whether covariant or contravariant. Even if 
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there is dimensional agreement, the numerical values may not be 
correct, because the base vectors are not unitary. 

Formulae have been derived that resolve this problem. If the 
physical components of a contravariant vector, A', are denoted by 
A(i), then the transformation law of physical components is 

A (j) = A ^ ( no sum on j) (2.29) 

In case of a covariant vector, first the associated contravariant 
vector has to be constructed and then equation (2.29) to be applied 

A 0) = A i ( no sum on j) (2.30) 

In addition, the transformation law of physical components 
between a curvilinear and a Cartesian system, is 

A curv. = /gjj A(i) cart < 110 SUm 011 J ) ( 2 - 31 ) 

By 


3. NAVIER-STOKES EQUATIONS 

As a starting point the momentum equations of Aris(ref.4), which 
are valid for a constant viscosity coefficient p, will be used 


au 1 

p -5T 


+ pu j u,j 


pf* + T* 


(3.1) 


where u'j is the covariant derivative of the velocity vector, f' is the 
vector of the external forces and T'ij is the covariant derivative of 
the stress tensor. The stress tensor for a Newtonian fluid is 
calculated from the following equations 


(3.2) 

(3.3) 
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T lj = (-p + Au,£) g' j + 2|ie y 
e ij = tfgV + gV")^ 


(3.4) 


These equations, e'i and e mn denote the deformation tensor in 
contravariant and covariant expression, respectively, and A is the 
bulk viscosity. 

After some manipulations, the stress tensor is found to be 
T ij = (-p + Au,£) g ,J + 1 1 ( g m ui + g* 1 u^) (3.5) 

For the calculation of the covariant derivatives, the following 
equations must be used 


i 



du 


dx 3 


f 1 * K 

7 + lj k> U 


(3.6) 


T* = J — £) + 1 * { ‘ k } ( 3 . 7 ) 

where J is the Jacobian and { j'|J is the second kind of the Christoffel 
symbol. 


If equation (3.6) is applied to the stress tensor, the following 
expression is found for its derivative 


T,f = g^(-p + AJ-2-(y)) + +g in ui» + 


dx? 


dx 


dx* 


( i . , Jm k kn j . 

< j U ’m + § U ’n ) 


(3.8) 


For the derivation of this equation, Ricci's lemma has been used, 
according to which the covariant derivatives of the metric tensors 
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vanish. Also, the rule that the covariant derivative of a scalar is 
equal to the conventional partial derivative, has been applied. 

In the absense of external forces, the momentum equations are 
given by 



j3u 

pu — + p{ 

ax’ 




(3.9) 


The continuity and the energy equations are easier to derive in 
generalized coordinates, because they can be expressed in 
divergence form. For example the continuity equation can be 
expressed as 


3a + V. (pc) = 0 (3.10) 

at 

and the energy equation can be expressed as 

3a + V. (ec) + V. (pc) = V.(T ij u. ) - V.q (3.11) 

at J 

where c is the velocity vector, e is the internal energy and q is the 
heat transfer. In the dissipation term, the stress tensor is 
multiplied by the velocity vector expressed in the covariant form 
because their product must be a contravariant vector. 

The heat transfer, in tensorian notation, may be written as 
q = - k VT = -kg' J — — (3.12) 

a^ 

Finally, the divergence of a vector a is expressed in tensorian 
notation as 


V. a = a, 1 . = J ( y a* ) 
1 dx J 


(3.13) 


16 


By substituting equations (3.12) and (3.13) in equations (3.10) 
and (3.11), the following expressions for the continuity and energy 
equations are obtained in the coordinate system x' = ( r\, £ ) 


i (A) + ^(py) + ±(p^) + 

at j J - T 


dr\ 


* 4 >-» 


(3.14) 


i.(£ )+ ± (( !lP)U) 
at J a^ J 


_a_ 

dr\ 


,e+P, 

J 


_a_ 


<(^)V) + — ((^)W) 

J dx 


k x 

8 jk u ) + 


+ -^-<Tg ij — ) (3.15) 

ax' J ax' 

In equation (3.15), the symbol x'i denotes the viscous stress tensor. 

In the next section, equations (3.9), (3.14), (3.15) will be used for 
deriving the boundary layer equations. 


4. BOUNDARY-LAYER EQUATIONS 

The boundary-layer equations will be derived in a system of 
coordinates, where the surface coordinates will be denoted by 
( x 1 , x 2 ) = ( rj ) , and the normal one by x 3 = t, . In the appendix it is 
shown that in this system, the metric terms g 13 , g 23 , gi 3 , g 2 3 
are equal to zero. Also, the Jacobian is simplified to the form 


j >/ £ 33^ S 1 1 S 22 ‘ §12)' 

For an order-of-magnitude analysis, the Navier-Stokes equations 
are expressed in a nondimensional form. So the freestream 
conditions, U^ p^ _ T„ a reference length L and the quantity 

p 00 U 2 oo (for energy and pressure) will be used. For these 
nondimensional variables, all the terms that include the viscosity 


17 


coefficients (n,X.) are multiplied by the inverse of the Reynolds 
number, based on the freestream conditions. 

If the nondimensional boundary-layer thickness is denoted by e, 
the normal velocity will be assumed to be W=0(e) and the Reynolds 
number 1/Re = 0(e 2 >- The metrics and the majority of their 
derivatives that appear in the Christoffel symbols, will be assumed 
to be of the order of 0(1). However, the normal derivatives of the 
metrics will be assumed to be 0(1) < gjj 3 < 0(1/8) , (i,j =1 or 2), so 
that a number of curvature-dependent terms will be included in the 
viscous part of the equations. Thus, it will be possible to obtain 
results of higher accuracy for bodies that have large curvature. 


4.1 Estimation of the Viscous Terms of the Boundary-Layer 
Equations 

In this section, the boundary-layer approximation will be applied 
to the viscous terms of the momentum equations. The viscous terms 
are 


x = 


Re 




0X* 


j < J <8 ‘Xi + S ,n U ’n » + * {/k }( ", 


kn j x 

+ g u >;) = 


= X + X, 


(4.1) 


Because of the presence of the Reynolds number, any term that is 
of the order of 1/e 2 , will be significant, so that the product of the 
term times the Reynolds number will be of the order of one. Thus, 
both, the partial derivative 3/3xi> and the covariant derivatives of 
equation (4.1) must be taken only along the normal direction 
(j,m,n = 3). Starting with the term x a ' for the momentum 
equation ( i=1) 


1 J ^ |L / 13 3 33 1 .. J d /U 


3 3 


1 V 

U n ) 


(4.2) 


The covariant derivative of the velocity component u 1 = U, is given 
by 
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U '3 = !t + { .‘3 ,U + , 2‘3 ,V + ( 3 1 3 |W “!t +I 3 1 1 |U+I 2 1 3 )V (4 ' 2) 
dQ dC, 

The Christoffel symbols that appear in equation (4.2) are given in 
the appendix . Each of them includes two of the normal derivatives 
of the following metric coefficients: g-n , g-j 2 .923. For example 

f 1 1 1 

M3' 2g ^22^33^11,3 * 2g g 3 3 g 2 1 %1 2 , 3 (4.3) 

Evidently, these derivatives depend on the curvature of the body 
along the and the t\- direction. For small values of the curvature, 
along one of these directions, the corresponding terms will be 
eliminated. 

For the t|- momentum equation ( i=2) 


2 ^ d /IX / 2 3 3, 33 2 .. J 3 , h 3 3 2. 

U '3 + S u ’3» = R?^r<fs “•:>> 


(4.4) 


where 


U '3 =|f + { 1 %}U + { 2 2 3 )V + {3%)W.|l + (, 2 3)U + { 2 2 , )V 
oC OS 


(4.5) 


The Tb' -component of equation (4.1) is estimated in a similar 
manner. For the ^-momentum equation(i=1) 


1 _LL r 1 i / 3 3 k k3 3. u , , 1 , 3 3 1,1.33 2. 

T b - Re ^3 k ^ u >3 + § u ’3 ) ~ p p ( {3 ilg u » 3 + (2 3^ u »3 ) (4-6) 


Re 


If the covariant derivatives are substituted in this equation, and 
terms including the products of the Christoffel symbols are ignored, 
as being smaller than 0(1/e 2 ), the following expression is found 



For the ri-momentum equation(i=2) 


t 2 = -W- / 2 H 33k 
T b Re i 3 k Mg ’3 


k3 3 . II , r 

+g 


3 3 


1 )g"u ,3 +( 2 2 3 )g"u,‘) (4.8) 


3 3 


Otherwise, after substituting the values of the covariant 
derivatives 


<- 


-1L ”, 
Re g 1 


2 \ aU 
3 l' 


a; 


JLg 33 ( 
Re 6 1 


2 

2 3 


} — 

ac 


(4.9) 


The Christoffel symbols that appear in equations (4.7) and (4.9) 
include the normal derivatives of the metrics along the surface 
coordinates. Thus, for bodies of small curvature, these terms may be 
totally eliminated. 

Finally, the term that depends on the bulk viscosity, X, has to be 
examined 


i i u a , 1T d V .. i ij a ,, u v w 

O = — g 1 — r ( AJ — - ( — )) =_g J _(XJ((_) +(_) +(_) )) 


Re 


3x k J 


Re 


ax’ 


J'§ J r\ J X 


(4.10) 


Since the dilatation terms are of the order of 1, no contribution 
from this term is expected. Furthermore, the term would become of 
the order of e for j=3. However, the metric g'i then becomes equal to 
zero, for i=1,2. Thus the bulk viscosity term will not appear in the 
boundary-layer equations. 

In summary, the viscous terms of the boundary layer equations 
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*■ ' + t* 33{ -'| < 4 ' 2 > 


4.2 Estimation of the Viscous Terms in the Energy Equation 


The dissipation term of the energy equation in a dimensional form 
is 


D = -^r(x lJ g jk u k ) = — 


Ox 1 


dx 


, , im j in i x k * 

•^<8 ». m + * "•„>Sjk u ) 


(4.13) 


If the term is nondimensionalized, the Reynolds number appears in 
the denominator. Thus again only the normal derivatives will 
contribute , i. e., i,m,n = 3. Next, the following approximate 
expression is found 


_ d ,33 j j3 3 . , 1 2 .v 

D= (p(g U , J 3 + g J U , 3 JCgjjU +g j2 U )) = 


d 331133 1 2 332 32 3 1 2 

— (^(g 33 U, 1 3 + g 1 V3)(g 11 U 1 + g 12 U 2 ) + ^(g 3 3 U, 3 2 + g 32 U, 3 3 )(g 21 U 1 + g 22 U 2 )) = 

<*» 


i-(Hg 3 V,(g l ,U + g 12 V) + ng 33 u.5(g ]2 U + g 22 V)) 


(4.14) 


The covariant derivatives of the velocity components are given by 
equations (4.2) and (4.5). Substituting them in the equation (4.14) 
and neglecting higher-order terms, it is found that 


^ a . 33, TT au au av V av u 

D =— (pg (gn U + g 12 V + gi 2 U + g 2 2 V ■))“ 

ac 11 ac 2 ac ac ac 
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±^iiL± 
ac 2g 3sac 


(U 2 ) + p- 


§22 a 

2g 3 3 ac 


2 S 17 9 

(V 2 ) + p-L2_ (UV)1 

g 33 ac 


(4.15) 


This is the final expression for the dissipation term. For its 
calculation, the metric relation : g 33 = g 33 has been used. 


4.3 Heat Transfer Term in Energy Equation 

The heat transfer term of the energy equation, in a nondimensional 
form, is given by 


H =-Vr gii 

ax 1 J 


aT a k uaT a k 2 j aT N a k 3j aT 

— ) = — ( -rS — ) + — (-r g — ) + — (-r g ~> 

dx‘ 3\ 1 3x J 3n 1 dx' a; J a x J 


(4.16) 


In the light of the approximation of the thin thermal layer , we 
retain only the normal-derivative term. If, the summation of the 
dummy index j is then applied , it is found that 


3 3 


a k 3 3 aT a ug • aa 

h = — (— g— ) = — (— ) 

ac J ac ac jpt(y-d ac 


(4.16a) 


where a is the local sonic velocity and the heat transfer coefficient 
has been substituted by means of the Reynolds analogy ( Pr = |ic p /k). 


4.4 Pressure and Other Terms in the Momentum Equations 


ine pressure 




i ij a p ii ap i2 ap i3 ap 

P =-g — =-g — -g — -g -r~ 


ax J 




an 




(4.17) 


The normal-derivative component is equal to zero , in both of the 
surface directions of the boundary layer ( i = 1 or 2 ) because of the 
zero value of the metric terms g 31 , g 23 . If , in addition, the 


contravariant metrics are replaced by the covariant ones, the 
following expressions are obtained 
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pi _ §3_3_£l_2 _ g 22 g 33 dp ^ ^ 

§ 3t 1 g at; 

p 2 _ g 3 3 g 12 3p _ g ll g 33 dp (4 19 ^ 

g a^ g ari 

Finally, there are the curvature terms, F', of the left-hand side 
(LHS) of the momentum equations 

F* = i.\) U j u k (4.20) 

For i=1 or 2, the curvature terms of the %- or of the ri-momentum 
equation are obtained: 

F 1 = { l \}U 2 +2{ 1 \}\]V + {^ 2 )y 2 (4.21a) 

F 2 = { 2 1 )U 2 +2{ 1 2 2 }UV+{ 2 2 2 }V 2 (4.21b) 

In the equations (4.21a) and (4.21b) the terms including the 
velocity component W have been eliminated, as they are of the order 
of e. 


4.5 Normal Pressure Gradient 


In the first-order boundary layer approximation, the 'n-momentum 
equation is commonly not considered , because a dimensional 
analysis shows that its terms are of the order of e . However in the 
present case, though the viscous and the convection terms are of 
this order, the curvature term is 0(1) < F 3 < 0(1/5) because this 
term includes the normal derivative of the metric coefficients 
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,3 , 3 , TT 2 . 3 


3 i,,2_ 1 3 3, i tt 2 ^12 TI „ , d %2 2 


F = { 1 }U +2{ 1 «}UV + L JV = - — g ( — — U +2 — — UV + — — V") (4.22) 

1 1 12 2 2 2 a; ac ac 


The corresponding pressure gradient is 


_ 3 3 3 3p 

P = ' g 

ac 


(4.23) 


So, in case of bodies of large curvature , a normal pressure gradient 
exists, given by: 


= + Uhl uv + l d ±l V 2 


ac 2 ac ac 


ac 


(4.24) 


4.6 Final Equations for Large Surface-Curvature 

If the boundary-layer approximations of the various terms 
.derived in this section, are substituted in the original Navier- 
Stokes equations of §3, the following higher-order boundary layer 
equations are found 

For the continuity equation 

7<f ) + ^ ( P7) + |-(Pt ) + ^ <p T ) = 0 (425) 

at J ac J ari J ac J 

For the momentum equations 
1. C- m omentum 


au TT au au w au _ tT 2 __ ¥TW .,2 

P “aT + PU TT + pV T" + pW 7T +K al pU + K a2P UV + K a3P V = 

01 ac ac 

K M - K ^ +~(t^-( — +K.JJ + K. ,V))+JLK. 8 ^. + JLK 0 ^- 


arj 


ac Re ac J S33 ac 


Re a8 


ac 


Re a * 


ac 


(4.26) 


2. ly-momentum 
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av av av av 2 2 

p - aT + pu TT + p v IT + pw ^7 + K b i pU + ^2 ?^ + k «p v = 
01 ac an di, 


K J2- - K ^ + ii(-^(-+K U + K V))+^K - + fK^ — 

b 4 -m-_ b 5 ar( Re a^ J g 3 3 dC, b 6 b 7 Re b 8 ;y Re b 9 




ac 


ac 


(4.27) 


and for the energy equation 


7 ( y ) + — ((> A> + — ((%W) = 
at J ac J dT\ 3 ac J 


1 ±(pK 2 &,K^) + 1 L^(-^ 


da 


Re 


ac 


ac 


ac 


e3 


ac 


Re 


ac Jgj 3 Pr(y -i) ac 


■) (4.28) 


These equations describe the most general case of a body highly 
curved in both the principal directions of the body axes. In this case, 
variation of the normal pressure gradient , equation (4.24), must be 
considered. 

The curvature coefficients appearing in equations (4.27) and 
(4.28) are defined by 


K , = 

al 



K a 3 ^2 2 ^ ’ 

K a 4 = 

g 33 g 12 

9 

g 

K a 5 = 

. g 22 g 33 
g 

K a 6 = 

la'il • 

K a 7 = ^2 3 ^ ’ 

xr 3 3 f 

^a 8 g ^3 

’ 

„ 3 3 

K a 9 g 

I 2 V 


K bl = 

{ 2 } 
l i it ’ 

cr 

N> 

II 

K> 

*S>3 = l'2 2 } ’ 

K b 4 = 

. g 33 g 1 2 

g 

K b 5 

_ g 3 3 g 1 1 

g 

K b 6 = 

l 3 1 > 

■ K b 7 =( 2 2 3 ) 

tr 33,2 , 

’ K b 8 ~ g I 3 1 1 

• K b 9 = 

3 3 r 

g ( 

2 I 

2 3 > 

K = 

Si 1 

k ,-iii 

, K = 

g 2 2 




el 

2g 3 3 

e 2 cr 

® 3 3 

’ e 3 

2g 3 3 





g 
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The Christoffel symbols appearing in these coefficients are given in 
the appendix. 


4.7 Final Equations for Moderate Surface Curvature 

In the case of bodies of moderate curvature, the terms including 
the derivative of the metrics in the normal direction may be 
omitted. Then the viscous terms of the momentum equations are 
considerably simplified. Only the first viscous term of the RHS has 
to be kept. Furthermore, in this case the metrics do not depend on 
the C -derivative . Equations (4.27) and (4.28) become 


P 


au 

at 


dU dU dU 2 2 

+ pu— +pv— ■ +pw— + K al pU 2 + K a2 pUV + K 3 pV 2 = 

dq drj dq 


K 


a 4 


dp _ K 

an’ ,5 di 


i a au 

+ — = — — (n — ) 

s 33 Re at at 


( 4 . 29 ) 


av av av av 2 2 

p ir +pu ^ +pv ^ +pw ^ +Kbipu +K b 2puv + K b3 pv = 


K ap ap 

K b4“ ■ ^ 


1 a av 

5 an + s 33 Re at ac ; 


( 4 . 30 ) 


In the first order of magnitude approximation it is sufficient to 
estimate the metrics and the Jacobian just on the surface (ref. 3), 
assuming that there is no significant variation across the boundary 
layer. The standard notation in this case is: g y= a ap , ( i,j,a,p = i , 2 ) , 

where the symbol "a" denotes surface quantities. Evidently, if this 
approximation is applied, the number of operations required for a 
numerical solution will be reduced considerably. 

The K factors that appear in equations (4.29), (4.30) are given by 
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K..-i,"*a + * g » (2 


28 35 


k 3 = ig 12 ^ + ts u < 2 — - — 
a3 2 an 2 


K = - g—U 

bl 2 8 35 


k =_ g ~_r^ 
b3 2 g ^ 


^12 

_ 

35 

an 

3 g, 2 

a s 22 

3 n 

" a* 

3 g, 2 


' ^ 

an 

^12 

^822 

an 

’ ^ 


). K ,2 

), K 


a 4 


li a 8ll 1 2 a &22 

g + g 

3n 


& 33&12 ^ & 22&33 

, iv c = 

g a5 g 


12 ^§1 1 22^§22 
+ g 


3 r| 

g 33&12 

g 


K b 5 = 


^1 1 ^3 3 
g 


If a fully-orthogonal coordinate system is examined, these 
relations are considerably simplified, because in this case the 
metrics gi 2 and g 12 become equal to zero. In addition, the covariant 
metric terms simply become the inverse of the covariant terms 
( g’i = 9ij ) and the Jacobian is simplified to: g=giig22933- If the 
metric terms are substituted by the hj-parameters(§2.1 ) and the 
physical-contravariant components are used 


h j = & ’ = AJ V^jj = AJh j (no sum on j) 

equations (4.29) and (4.30) simplify to 


3 U 


u au 


p-rr- + p— + P 

3 t h i 35 


v au 

h 2 arj 

au 


+ P - 


w au 
l 3 ac” 


+ pUVK t - pV K 2 = 


1 a P 1 a , 

T + “= (P — ) 

h i ^ h?Re3C ^ 


(4.31) 


av 

at 


u av v av w av 2 

+ + p ““ + p uvk 2 - p u K i = 


+ P h i a^ 


2 aii 

1 a P 1 a , av 

‘ FT + ~1 — — (PTT) 
h 2 arj h^RedC a C 


*3 K 


(4.32) 


where K-j, K 2 denote the geodesic curvatures of the surface lines 


27 


K 1 = T—r 1 , K 2 = h-h ^ 

h 1 h 2 3rj h l h 2 a$ 

In these well-known equations, the basic contravariant symbols 
(U, V, W) have been used for denoting the physical velocity 
components. 

In the appendix, the first-order boundary layer equations are given 
in the case where the normal to the surface is denoted by x 2 = r\. 


5. CONCLUSION 

A set of higher-order boundary-layer equations have been derived 
for a nearly-general curvilinear coordinate system. In this system, 
the surface coordinates can be nonorthogonal, whereas the third axis 
is restricted to be normal to the surface. For the derivation of the 
equations, an order of magnitude analysis has been performed in the 
Navier-Stokes equations, written in a fully nonorthogonal coordinate 
system using the tensor notation. Assuming that the normal 
derivatives of the metrics are: 0(1) < g ( j 3 < 0(1/8) , (i,j =1 or 2), a 
number of curvature-dependent terms appear in the viscous part of 
the equations. Thus, it will be possible to obtain results of higher 
accuracy in the case of bodies that have large curvature. 
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APPENDIX 


Metric Tensors 

By application of the definition it is found that the elements of 
the covariant metric tensor are given by: 


g xl = x\ + y\ + z\ , + + . «i 3 -Vc + y *W; 


% $ 

2 2 2 

g = x + y + z 

1 1 T | T | Tj 

2 2 2 
§33 = X ; + 


\ T | % 11 

g,, = x x +y y +z z 
23 ii c ii c ii ; 


(a. 1) 


The elements of the contravariant metric tensor can be estimated 
similarly by simply interchanging in equations (a.1) the Greek 
letters ( % £ ) by the Latin letters (x, y ,z). Alternatively they are 

given as the inverse elements of the covariant metric tensor 
(equation 2.13) 

2 

11 _ §22 §33 " §23 


12 §23 § 3 1 ' § 21 §33 

■ — i — 




22 

s 


§11 §33 ~ §13 

s 


3 3 

g 


2 

§1 1 §2 2 ' §12 

g 


23 


g 


§21 §31 ~ § 1 1 §32 

g 


(a.2) 


Note that both of the metric tensors are symmetric. 

Some of the elements of the metric tensors become equal to zero 
in the case of the partially orthogonal system assumed in this paper, 
because the following relations are valid between the base vectors 
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e 2 e 3 - % ■ x ; + % + % ^ = 8j3 = 0 (a3) 

Thus, the metric tensor and the Jacobian are simplified 
considerably 



By considering the relation that connects the covariant and the 
contravariant metric tensors it is easy to show that the 
contravariant metric terms g 13 and g 23 are also equal to zero. 

Christoffel Symbols 

The Christoffel symbols are given here only for the particular 
semiorthogonal coordinate system examined in § 4 


1 , 1 11 a §ll . 1 , 13„ ^13 13 a §ll x f 1 , _ 1 11 3g ll A 1 13 


2 * s an 2 - ^ 


) . U 2 ) = js 


+ T § 


13 


f 1 l--(g n ^-+g 13 ^) {' } = — g 11 ^11 + — g 13 ^ §33 

ll3} 2 (g ac +g ^ } ’ 2 g ^ 2 S 


an 


, 1 , 1 „13 ^33 , 1 ,JU 9g 13 J1 a %3 , ,3 , _ 1 13 d &U 1 33 ^13 

3 3 } = "2 g l^ + 2 (g 2 — g 17 + 2 g a. 


ac as 


, 3 , 1 13 ^11 1 , 33 ,, 5 gi 3 33 ^§11 . , 3 . _ 1 13 ^13 1 33 ^33 

i i = T g l^ + 2 (g 2 *" g ir )’t 23 }-2 g ^ + 2 g 


a$ ac 
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3 , 1 _33 33 1 13 9 ^®13 13 ^33 .3 ,1 13 ^§11 1 33 ^§33' 

ac 2 ac ^ 1 3 2 ac 2 as 


^ 3 3^ ~ ” 8 


f 2 1 1 22 ^Sll ,2,1 22 ^§22 f 2 l _ 1 22 ^8 13 

‘i >’ 2 g an ' 12 2 a$ ' 13 2 an 


1 22 ^§22 


{ 2 v=i g ~-^- ■ * t 3 v=-Tg 


1 _22 ^22 
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2 , 1 _ _22 ^833 

2 ' arj 


The Boundary Layer Equations in an Alternative System 

If the normal direction is designated by x 2 ■ ti and the crossflow 
direction by x 3 = £, the boundary layer equations in the case of the 
first order of approximation become 

for the ^-momentum 


P 


au 

at 


dU dU dU 2 2 

+ pU +pV -_ + pw— + K lP U 2 + K a2 pUW + K 3 pW 2 = 

dq dn dL, 


K <>? . L_ JLf m — t 

a4 ac a5 a^ 8 2 2 Re an ail 


(a.5) 


for the ^-momentum 


aw aw , aw aw tt2 Tjr ttw Tjr „ r2 

P 1T + pu 7T +pv TT +pw 77 +K biP u +K b2 puw + K b3 pw = 

dq dn oC, 


K b4f-Kb5^7 + 1 

aq 


_a_ aw 
a; 822 Re an an 


(a.6) 


In this case the K values are given by 
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